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Abstract 

Backgrounds are pervasive in almost every application of 
general relativity. Here we consider the Lagrangian for- 
mulation of general relativity for large perturbations with 
respect to a curved background spacetime. We show that 
Ncether's theorem combined with Belinfante's "symmetriza- 
tion" method applied to the group of displacements pro- 
vide a conserved vector, a "superpotential" and a energy- 
momentum that are independent of any divergence added 
to the Hilbert Lagrangian of the perturbations. The energy- 
momentum is symmetrical and divergenceless only on back- 
grounds that are Einstein spaces in the sense of A. Z. Petrov. 



Here we use well tested methods in classical field theory to construct 
a conserved vector density with respect to arbitrary backgrounds. 
Backgrounds are pervasive in almost every physical application of 
general relativity: from gravitational radiation [1] to testing the 
theory in the solar system at higher orders of approximation in the 



PPN formalism [2] ; from stability theory of de Sitter or anti-de Sit- 
ter spacetimes [3] to stability of black holes [4] . Relativistic cosmol- 
ogy is studied on a Friedmann- Robertson- Walker background [5]. 
Isolated sources are analyzed on asymptotically flat backgrounds. 
Thus, backgrounds play an important role in practically all ap- 
plications of general relativity. It is therefore interesting to have 
satisfactory differential conservation laws on curved as well as on 
flat backgrounds. 

There are essentially two methods to obtain conserved vectors 
in general relativity. One method consists in rewriting Einstein's 
equations for the perturbations keeping on the left hand side terms 
linear in second order derivatives of the perturbed gravitational field 
components ([6] [7] on a flat background; [3] on a Petrov space [8] 
background). Einstein's equations have been obtained in that form 
directly from a variational principle [9] [10] [11] [12]. The right hand 
side of the equations is a symmetric "energy-momentum tensor" say 
*T^ U and for any Killing vector of the background £„ there exist a 
conserved vector density ^/—g{*T^ v )^ v . 




There are problems with that method. First, as pointed out 
by Boulware and Deser [13] and by Popova and Petrov [14], the 
perturbations of the gravitational field can be represented with the 
metric, the inverse metric, the metric density and so on. For each 
representation the conserved vector density will be different. Sec- 
ond conserved vectors have always been obtained for Killing vectors 
only and only on A.Z. Petrov spaces as backgrounds. 

A second way of finding conservation laws consists in applying 
Ncether's method to a Lagrangian of the gravitational field (see for 
instance [15]). This leads to a "canonical" Ncether conserved vec- 
tor, the divergence of a superpotential and to a canonical energy- 
momentum tensor. The method gives conserved vectors on any 
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background with any vector field defining a one parameter displace- 
ment [16], not only for Killing vectors of the background. This 
great multiplicity of conservation laws in general relativity is re- 
lated to the relabeling of spacetime points and is not without anal- 
ogy with circulation conservation in fluid dynamics. The relabeling 
of points in barotropic flows is associated with Ncether conserved 
vectors which in comoving coordinates are known (not too well) as 
conservation of "potential vorticities" [17] [18]. The conservation of 
potential vorticities is the local expression of the more familiar non 
local Kelvin conservation law of circulation. 

But there are problems with this method also. First the La- 
grangian density is not unique. A divergence can and must be 
added to the Hilbert Lagrangian because the Hilbert Lagrangian 
leads to Komar's [19] conservation law which gives the wrong mass 
to angular momentum ratio in the weak field limit (the "anomalous 
factor 2" [20]) and does not give [21] the Bondi [22] mass. Various 
divergences have been added for different reasons [23] [24] [25] . They 
lead to different conserved vectors. Second, the canonical energy 
momentum tensor is not symmetrical nor in general divergenceless. 
On a flat background the energy-momentum is divergenceless but 
it is not symmetrical and does not provide a conserved angular mo- 
mentum and when it does, the angular momentum does not include 
the helicity of the field. Thus this second method is not satisfactory 
even in the weak field limit. 

Here we use Ncether's method, the Hilbert Lagrangian density 
but also Belinfante's [26] modification. In classical field theory, 
Belinfante's correction gives a symmetric field energy-momentum 
tensor which ensures conservation of angular momentum, helicity 
included. As we shall see, the Belinfante correction has the great 
advantage to provide a conserved vector and a superpotential that 
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are independent of any divergence added the the Hilbert Lagrangian 
and does not depend on a particular representation of the gravita- 
tional perturbation. The result looks like a peculiar blend of the 
two methods with none of their defects. The answer is unique. 

Let us start with the Lagrangian for gravitational perturbations, 
the matter Lagrangian playing no role in our considerations:^ 

^ = -^i k - t ) = iy -^- « 

R = y/—gR is the scalar curvature density of a spacetime with 
a metric g^ u , R is that of a background with a metric g^, both 
metrics have signature —2, and k = %itG/c A . A hat "~" always 
means multiplication by y/—g not by y/—g, thus R = <>/^gR is 
different from R = y/—gR. 

We now apply Ncether's theorem to £', not CJ G . For this we first 
calculate the Lie derivative £^C of C for an arbitrary displacement 
field it is equal to (C'^J and is of the form 

£^ = j- K & V£ t^ - j-^ W - §r£^) = (c'e) 

(2) 

where G^ v is Einstein's tensor, the left hand side of his equations. 
We then use the contracted Bianchi identities D V G^ U = and Ein- 
stein's equations G^ = nTj^ and obtain a conserved vector density 
i ' M which looks as follows 

d,i = o, t = tx~y k {g pa £^ - r p £^r; a ) (3) 

We then redo the same calculations with the Lie derivative £%£ 

of C' which is equal to d^C'^) and we obtain a conserved vector 

1 Symbols may seem unnecessarily complicated for this short paper but 
they are the same as in the full paper where they have some justification. 
Comparison will thereby be simpler. 
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density i that satisfies exactly the same equality (3) with bars 
over every symbol except C, p . 

The conserved Ncether vector density J' M associated with CJ G = 
£' — C' is equal to the difference t'^ — t 'v. This difference is a linear 
homogeneous expression in in ~g -covariant derivatives D p ^ a 
(with £ CT = g^^) and in derivatives of -D(p^o-) which we denote by 
a special symbol z pa ; thus J' M can be written in the following form 

where l pa = g pa — g pa ; indices are displaced with g^, never with 
g pv . In Eq. (4) §' p is the relative energy momentum tensor density 

Q,tx = f» _f» + a f^id comp onent. (5) 

The field component of the tensor is complicated as may be guessed 
from Eqs. (3) and, say, (3); however, its explicit form will not be 
needed and there is no point in writing it here in detail. The a' wu 
term is more important. Its antisymmetric part a' p ^ plays the 
role of a relative helicity in linearized quantum gravity [27] and is 
similar to the helicity in electromagnetism 



2ko»*° = g^r^ix + g uX r a Kx - zg up r x Kx (6) 

where the tensor = Y^ v — is the difference of the Christoffel 
symbols. 

It is well known [29] that the conserved vector I' p is equal to 
a divergence of an anti-symmetric tensor density, a superpotential, 
K pv which for good reasons we call the "relative" Komar [19] su- 
perpotential, relative to the background 

I' p = + b'w°D p i c + fj p = d u K p \ 

K pv = -K up = ± (gP^Dp^ - gPfrDpfl) . ^ ' 
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contains g M ^-covariant derivatives D p as well as g^-covariant 
derivatives D p . 

We now apply Belinfante's [26] procedure. The method has been 
used in general relativity by Papapetrou [30] on a flat background 
with Killing vectors of rotation in Minkowski coordinates. It is here 
applied to curved ones with arbitrary £'s and in arbitrary coordi- 
nates. It works as follows. Replace the conserved vector density 
by the following new (also divergenceless) one X M 

In the new vector, there are no antisymmetric derivatives of £ any- 
more, the S"-addition cancels precisely the helicity-term. The S'- 
addition also modifies the energy- momentum tensor density 9 f Jf, the 
ff- vector density and the superpotential K^ u . J M has the following 
form 

The new energy-momentum tensor density TJ 1 and the new Z^ are 
related to 9'Jf and f]^ as follows: 

= 8* + DpS'^, Z" = ff+ + S'^ a )z pa (10) 

while the Komar relative superpotential is replaced by a new su- 
perpotential 

= + S'^ p ^. (11) 

Notice that Z^ like f]^ is zero if £ M is a Killing vector ^ of the 
background. 

One crucial point is now this: let's add a divergence to R — R in 
the Lagrangian density (1) say d^. This has the effect to produce 
another Ncether conserved vector density 1^ ^ 1'^. Indeed, the Lie 
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derivative of the divergence, £^{d IJi k lx ) = d^(£^) = du^d^k^). 
Since £^k fJ ' contains at most first order derivatives of £ M the con- 
served I 11 will have a modified £ M factor (9£ ^ 9'Jf) and a modified 
D p ^ a factor (a Mpcr ^ d' wa \ Of course the superpotential is also 
changed and it is easy to find that K^ v is replaced by 

= + (12) 

With a different <7 MP<:r there is also another S fiup as can be seen from 
Eq. (8) and it is equally easy to find how that S is related to S': 

S^p = S'^ p --^H u] . (13) 

K 

Thus, look at Eq. (11), X MZ/ does not depend on adding a divergence 
to the Hilbert Lagrangian.f\ This is also true for TJ 1 and Z^ as 
implied by Eq. (9). 

The explicit structure of 7~f is not important here; it can be 
derived from Eq. (15) given below . What is important however 
are the properties of TJ 1 which we most easily obtain by Rosenfeld's 
[32] method. The modified conservation law <9 M Z M = is linear in 
£ M with derivatives up to order three that come from d^Z^ as can 
be seen from Eqs. (9), (10) and (4). dpi** = may thus be written 
in the form 

dfi» = fa" + P^D.C + P P \D (P «)C + P^uDqv*)?' = 0. (14) 

This equation holds for arbitrary Therefore all the /3's must 

be equal to zero. These are the "Rosenfeld identities". The most 

interesting identities for now are those involving TJ 1 which, we can 

2 Bak, Cangemi and Jackiw [31] made already the interesting remark that 
Belinfante's modification of the Noether currents obtained from Hilbert's or 
Einstein's Lagrangians lead to the same symmetric and divergenceless energy- 
momentum tensor relative to a flat background in Minkowski coordinates. 
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see from Eq. (9), are /3^ v = or rather = ^ p g pv = and 
Pv = 0: 

piiv = fnv + 7^ hji(.pny + &{in>)n - a'(Hp} - ^ X R U X = 0, 
and 

A, = - ^D u R p(7 = 0. (16) 

The new energy momentum tensor density has the following form 

f ^ = (t^p - TiA + f pu + - [\l prj R p(T g pv + [ A Kft^ . (17) 

It contains three types of terms, a symmetric matter energy-mo- 
mentum of the perturbations, a symmetric field energy-momentum 
tensor| T pv = r up and two non-derivative coupling terms of the 
metric density perturbation to the background Ricci tensor, the 
last one only being anti-symmetric in fiu. Therefore: 

(i) Eq. (15) shows that f pv = f vp if and only if l x ^R^ = 0, 
i.e. if R^ v = —A.g with A necessarily constant. Thus the energy- 
momentum tensor is symmetrical only if the background belongs 
to the class of Einstein spaces in the sense of A.Z. Petrov [8]; these 
are not the "Einstein spacetimes" of Friedmann- Robertson- Walker 
cosmologies. However de Sitter and anti-de Sitter spacetimes are 
Einstein spaces. 

(ii) For those particular backgrounds, Eq. (16) shows that T pu 
is also divergenceless: D u T pv = 0. There are thus no divergence- 
less and symmetric field energy-momentum tensors of perturbations 
except on Einstein spaces. 

(iii) If the background is an Einstein space and £ v one of its 

Killing vectors, the conserved vector density has the simplicity of a 

classical expression J p = T pl/ ^ v = d u J pv . 

3 This field energy-momentum tensor is a horrifying quadratic homogeneous 
expression in F a , Dj pa and D^IP". 
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(iv) Notice that if the background is not an Einstein space and 
£ M is not a Killing vector of the background there are still plenty of 
conserved vectors as can be seen in Eq. (9). 

Having considered the properties of we now turn our at- 
tention to the superpotential which has a rather simple form on all 
backgrounds and with every vector 

This superpotential generalizes a number of well known particular 

cases: 

- On flat backgrounds with Killing vectors of rotations and in 
Minkowski coordinates, T^ u is the superpotential found by Papa- 
petrou [30]. V^ up has been occasionally refered too (see for in- 
stance [33]) as Papapetrou's superpotential. This is strictly true 
with Killing vectors of translations only (and on a flat background 
in Minkowski coordinates). 

- The tensor density V^ vp is the same expression as Weinberg's 
[6] Q^ p and Misner, Thorne and Wheeler's [7] d a H^ vp though in 
those expressions is the linear approximation of the inverse of 
hpu 9pu 9/j-u- 

- P^p is identical with the linear approximation on a flat back- 
ground of Freud's [34] superpotential, and of Landau and Lifshitz's 
[15] superpotential. 

- The complete superpotential (18) with Killing vectors, j^, is 
similar to that of Abbott and Deser [3]. To obtain their superpo- 
tential replace l^ by -yf^g times = H pa g p ^g au with H pu = 
h pu — ^9 P u9 pct ^p<t- I n the linear approximation \p v = —y^gHP u ; 
the two superpotentials are thus (also) equal to the lowest order 
in h^ u , but not to higher orders. Both superpotentials give the 
total energy and angular momentum for stationary spacetimes at 
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spatial infinity. But, the superpotential (18) alone gives the four 
momentum of Bondi [22] and Sachs [35] 's radiating spacetimes at 
null infinity. 

Eq. (11) can be looked at as a "corrected" Komar superpo- 
tential. There are other corrected Komar superpotentials in the 
literature [36] [37]. These modified Komar superpotentials have the 
anomalous factor of 2 mentioned above and have also other unsat- 
isfactory features [38]. 

Differential conservation laws on a curved background are useful 
in relativisitic cosmology and have indeed already been used. For 
instance, Friedmann-Robertson- Walker spacetimes admit a time 
translation conformal Killing vector. The corresponding conserved 
current helps solving Einstein's equations with scalar perturbations 
and topological defects in the limit of long wavelength [39] [40]. An- 
other instance is the "integral constraint vectors" which were used 
by Traschen and Eardley [41] to analyze measurable effects of the 
cosmic background radiation due to spatially localized perturba- 
tions. Those Traschen [42] vectors can be shown to be linear com- 
binations with cosmic-time dependent coefficients of conserved vec- 
tors associated with the conformal Killing vectors of "accelerations" 
described in Fulton, Rohrlich and Witten [43]. 

We believe that finite volume integrals of the conserved vector 
densities, which are equal to closed surface integrals of the superpo- 
tential may be useful in numerical calculations for the same reason 
that the relativistic virial theorem of Gourgoulhon and Bonazzola 
[44] is useful to check numerical integrations of relativistic neutron 
stars [45]. 

Detailed calculations not included in this letter will appear in a 
full paper presumably in Class. Quantum Grav.. 
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